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A NUMERICAL METHOD FOR 
UNSTEADY AERODYNAMICS VIA ACOUSTICS 

STEVE HODGE* 

Abstract. Formal solutions to the wave equation may be conveniently described within the 
framework of generalized function theory. Here generalized function theory gives a formulation and 
formal solution of a wave equation describing oscillation of a flat plate from which a numerical method 
may be derived. 


Summary. Wave equations describe vibrations and spatial perturbations (“waves” ) 
of physical terms away from certain ambient terms. In acoustics vibration terms are 
generally described as “sound” though they are of a wider range than audible sound. 

In general, think of wave equations as describing any moderately small variation in 
space or time. Acoustics and aeronautics parallel in that the governing equations of 
linearized compressible aerodynamics and acoustics are the same. The linear acoustic 
wave equation for the velocity potential <j>{x,t) is 


(i) 


dt 2 c 2 dx 2 " 


where c is the “speed of sound.” The nature and scope of the infinite possible solu- 
tions to (1) vary enormously. Usually a particular physical situation or other specifies 
“initial conditions” or “boundary conditions” which narrow the admissible solutions 
down to one. However, there is an enormous gulf between knowing that a unique 
solution is there and writing it down with a particular equation or function from 
which to calculate. There are many ingenious techniques for bridging this gulf. In 
particular, sound generated from physical situations which involve moving objects 
may be described by wave equations with “generalized derivatives” and have formal 
solutions which involve integrals of generalized functions. Generalized derivatives are 
a useful formal method roughly equivalent to a “weak” formulation [1], but with the 
benefit of derivative formalism in which the work of integration by parts is automati- 
cally incorporated. The generalized functions and derivatives are useful in describing 
discontinuous phenomena such as the intrusion of a wing in free space or a “shock” 
in flow quantities on the surface of the wing [2]. 

In thin airfoil theory the boundary conditions are specified on the mean chord 
surface [3]. In terms of generalized derivatives and linearized momentum equation 
p = —p Q 2& y the wave equation becomes 

(2) ff - ^0 = | \po v n | V f me | A(/ mc )] + V • [A pVK6(K)) 


where f mc is the mean chord surface, K is the mean chord surface and wake, v n is the 
normal velocity due to thickness and the bars over the derivatives denote generalized 
derivatives [4]. This equation is similar to the Ffowcs Williams-Hawkings (FW-K) 


(3) W - = Ft + 9 • b*W)l 
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which is a formal restatement of the equations of continuity and momentum [5]. Since 
pressure is continuous across the wake, this may he rewritten 

(/|) = ji \P0Vn\Vfme\MImc)] + V • [ApV/ me «(/ me )] . 

Equation (3) lias formal solution 


(5)ljrp(x,t) 


\d_ f f PQCVn ± Apcos(0) \ r / Apcos(6) \ 

cdtJ Jmc=l) \ r|l — M r \ ) T . ' J jmc=0 \r 2 \l- M r \) T . 


dS 


Here r* is the emission time of a signal received by the observer at time t. Integration 
of (5) with respect to i and the identity p = —pojj- t <i> leaves 


4W(«,0 = -'/ dSdf. 

Cjj mc =„\ r| 1 — M r \ ) T . \r*\] - M r \J T ., 


( 6 ) 

Here t*' is the emission time of the signal received at time t' . 

Generally (6) is separated into two terms. We are interested in the so called 
“lifting” terms (corresponding to aerodynamic lift), namely: 


(nyp„4>(x,t) 


-H. 


/ A pcos(0) \ 

=o V r 1 1 — Af r |/ r . 


dS + 


f j 

J-oo JUo - 0 


/ A pcos(0) \ 

\r 2 |l — A/ r |/ T ., 


dSdt'. 


Refer figure 1 for the following. We assume that the mean chord surface is a flat 
plate moving with uniform velocity in direction shown. Assume that the pressure 
is piecewise of the form A pe ,r< over individual surface panels of the plate. This 
corresponds to an unsteady aerodynamic loading. On individual panels in the 
coordinate system shown in figure 1 the individual integrals restated in terms of 
normal velocity reduce to 


\xv n (x,l) = 1/c— — 
ux:\ 

( 8 ) 


f / Apcos(O) \ 

],•„ V\\-M r \J T . 


dS + 


a 

dx:\ 


[‘ f ( A pcos(0) \ 

W\l- M r \) T ., 


dSdt' 


where v n is the velocity in the normal direction and M r is the Mach number in the 
r-dircction. Since wc are interested in aerodynamic lifting, the observer point will 
be on the lifting surface and consequently (8) will be singular on panel on which 
the observer lies. Although divergent, a physically correct solution may be extracted 
when the integral is considered in the Cauchy principle value sense with corresponding 
“Hadamard finite part” [()]. After a good choice of coordinate systems, reducing to 
a limit definition, using a principle value, and alot of summertime work, the integrals 
remarkably reduce to lim/,_>o+ + h) where 


(9) /i 

and 


= h/cc iVi I 
Jr i 


r j - — iTr / c 


f . -1 r y-j(r) 1 

r m [v/p^J 


IVr 2 - /» 2 



(the endpoints r, are the range of r over pane) l\j and <|uanlitics in the 2nd integral y, 
ami r depend on t' ). We arc now ready to consider the global effect of simultaneous 
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observers in each panel. This leaves a system of algebraic equations which may be 
solved for unknowns A p. At present a numerical method for this system is being 
formed with the aid of Mathematica. 
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